Elastic moduli renormalization in self 
interacting stretchable polyelectrolytes 

Rudi Podgornikf* Per Lyngs Hansen and V. Adrian Parsegian 

Laboratory of Structural Biology 
National Institute of Child Health and Human Development 
National Institutes of Health, Bethesda, MD 20892-5626 

and 

f Department of Physics, Faculty of Mathematics and Physics 
University of Ljubljana, SI-1000 Ljubljana, Slovenia 
and Department of Theoretical Physics 
J. Stefan Institute, SI-1000 Ljubljana, Slovenia 

February 1, 2008 



Abstract 

We study the effect of intersegment interactions on the effective bending and 
stretching moduli of a semiflexible polymer chain with a finite stretching modulus. 
For an interaction potential of a screened Debye-Huckel type renormalization of the 
stretching modulus is derived on the same level of approximation as the celebrated 
Odijk-Skolnick-Fixman result for the bending modulus. The presence of mesoscopic 
intersegment interaction potentials couples the bending and stretching moduli in 
a manner different from that predicted by the macroscopic elasticity theory. We 
advocate a fundamental change in the perspective regarding the dependence of 
elastic moduli of a flexible polyelectrolyte on the ionic conditions: stretchability. Not 
only are the persistence length as well as the stretching modulus dependent on the 
salt conditions in the solution, they are fundamentally coupled via the mesoscopic 
intersegment interaction potential. The theory presented here compares favorably 
with recent experiments on DNA bending and stretching. 



*To whom correspondence should be addressed. E-mail address: rudi@helix.nih.gov. 







1 Introduction 



Mechanical properties as characterised by bending, stretching and twisting and their 
respective elastic moduli, more than anything else determine the supercoiling [Q and 
packing properties of DNA ||. The effect of intersegment interactions along a semiflexible 
polymer chain on its bending modulus, especially in the context of DNA, has been studied 
for quite a while (see Ref. ||] and references therein). The major result of these studies 
is the celebrated Odijk-Skolnick-Fixman (OSF) formula [|J that connects the value of the 
persistence length with the parameters of the interaction potential between the segments 
along the polymer chain usually assumed to be of the Debye-Hiickel form. 

In the case of a stretchable chain there exists no theory that would connect the inter- 
segment interactions and the stretching elastic modulus. That this theory is very much 
needed is shown by recent experiments on single molecule DNA stretching ||. They 
point to the conclusion that the measured renormalization of the stretching modulus in 
the presence of solutions with different ionic strengths can not be rationalized on the 
basis of simple elastic theory arguments || according to which the renormalization of 
the bending and stretching moduli should be proportional. Experiments on the contrary 
suggest that while the presence of electrostatic repulsion between DNA segments tends 
to stiffent up the chain it also makes it a lot more stretchable. This would suggest that 
DNA might not be describable by macroscopic elasticity theory at all |]5[. 

The physical basis of the bending modulus renormalization is quite simple and well 
understood, see Fig. ??. The change in the persistence length is due to the fact that the 
effective spacing between the segments gets smaller (L' < L) after the polymer is locally 
bent. Because the intersegment interactions are assumed to be repulsive the interactions 
oppose bending and thus give rise to a higher bending modulus. The physics of the effect 
of the intersegment interactions on the stretching modulus, being proportional to the 
second derivative of the interaction energy as a function of the intersegment coordinate, 
is quite different. If the chain is allowed to stretch locally the segment length becomes 
bigger (a' > a), the interaction energy becomes less steep and its second derivative thus 
becomes smaller. Therefore the more the segments are further appart the less the second 
derivative of the interaction energy is going to be, leading to a smaller stretching modulus. 
In what folllows we will formalize and extend this simple physical picture. 

One should note here that part of this effect already transpires through the work of 
Ha and Thirumalai 0. Though they deal with a nominally unstretchable chain, the 
unstretchability coinstraint is implemented globally through an appropriate Lagrange 
multiplier. In the presence of intersegment interactions this Lagrange multiplier would be 
renormalized by the interaction. Our calculation builts on and adds to this change in the 
Lagrange multiplier by explicitely introducing a stretching part of the elastic energy. 

In this article we present a straightforward generalization of the OSF arguments to 
include the effect of the intersegment interactions on the stretching modulus of the chain 
as well. We derive both the bending modulus or equivlanetly the persistence length renor- 
malization as well as the stretching modulus renormalization concurrently on the same 
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level of approximation based on the recent implementation of the 1/d expansion tech- 
nique into the self-interacting semiflexible polymer theory ||. We show that for finite 
range intersegment interaction potentials the bending and stretching renormalization be- 
come coupled. Comparison with recent experiments seems to bear out our line of thinking 
quite strongly. 

We are convinced that our calculations dispell any doubts as to whether DNA conforms 
to macroscopic elasticity theory ||. It does, if one takes the long range part of the 
interaction potential between the segments of the polymer chain consistently into account. 

The outline of the paper is as follows: First we introduce a mesoscopic model of the 
self-interacting polymer chain with bending and stretching elasticity included. We briefly 
describe the 1/d expansion method that we use to evaluate the partition function of the 
model. We then explicitely obtain a mean field solution of the model and the contribution 
of the fluctuations to the equation of state of the polymer chain that connect the elongation 
of the chain with the stretching force acting on it. We finally derive the renormalized 
elastic moduli of the chain and compare the results with available experiments. 



2 Model 

We start by formulating an elastic mesoscopic Hamiltonian for a self interacting chain 
described in the highly stretched, small deformation limit in the Monge-like parameteri- 
zation (see Fig. ??) as r(s) = (z,p(z)), where s is the arclength along the chain. For a 
one-dimensional solid |J which is an adequate representation of a flexible polymer chain 
on small lengthscales, the deformation tensor has only one non-zero component 

du z {z) l ( dp{z) \ 2 

where u z (z) is the internal phonon-like field describing the stretching of the chain. The 
bending field p(z) is in the direction perpendicular to the local tangent of the chain, thus 
perpendicular to the z axis. This result can be derived straightforwardly from the form 
of the line element along the chain: ds 2 = (dz + du z (z)) 2 + (dp(z)) 2 0. 

The total mesoscopic energy of a semi-flexible self interacting chain under external ten- 
sion contains three terms: the stretching elastic energy term, the bending elastic energy 
term, the stretching force term and the interaction term. It can be written straightfor- 
wardly in the form 

- // ^^-dz + \f J V(\r(z) - v{z')\)dzdz\ (2) 

where V(|r(z) — r(V)|) is the interaction potential between two segments of unit length 
separated by |r(z) — r(Y)| 2 = (z — z' + u z {z) — u z {z')) 2 + (p(z) — p{z')) 2 . f is the external 
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force stretching the chain in direction z, A is the stretching (Lame) modulus and Kc is 
the bending modulus related to the persistence length l P as K c = kT l P . This mesoscopic 
energy presents a generalization of the existing models of a stretched elastic chain, Refs. 
and [TJ. 



3 The 1/d Expansion Method 

The non-local nature of the intersegment interaction, dependent on both z and z', pre- 
cludes an explicit evaluation of the partition function of the model with a mesoscopic 
Hamiltonian Eq. |||. We thus have to resort to some approximation scheme that will make 
the evaluation of equilibrium properties of this model tractable. 

At this point we introduce the 1/d expansion method to obtain an approximate but 
nevertheless explicit form of the partition function. The basis of this method is the 
introduction of two auxiliary fields: B(z, z') = (r(z) — r(V)) 2 and its Lagrange multiplier 
g(z,z') (for details see ||) that will help transform non-local intersegment interactions 
along the chain into local energy terms. With these variables, and limiting ourselves to at 
most quadratic order in all the variables, the chain mesoscopic Hamiltonian can be cast 
into the form 

n = ^/(^) 2 ^ + |^/(^) 2 ^-/ML)-^(o)) + 

+ \J J dzdz'V{B{z,z')) + 

+ 1 J J dzdz' g(z, z') (B(z, z') -(z-z' + u z {z) - u z {z')f + (p(z) - p(z')f) , 

(3) 

where we indicated explicitely the dependence of the interaction potential on the auxiliary 
field B(z,z') as V(B(z, z')) = V(^B(z, z')). This dependence is introduced via the 
Lagrange multipier in the last line of the above equation Eq|3] through the constraint 
B(z,z') = (r(z) -r{z')f. 

The rationale for this change of variables is that the dimensions of the fields p(z) 
and u z (z), that can be integrated over explicitely and exactly, are assumed to be much 
larger than the dimensions of the auxiliary fields B(z,z r ) and g(z,z f ). This allows the 
contribution of the auxiliary fields to the partition function to be evaluated on the saddle- 
point level. This approach can be shown to be asymptotically exact if the dimension of 
the embedding space for the polymer chain, i.e. the dimension of the r(s) vector, tends 
to infinity. If this is not the case, as in deed it is not for our three dimensional case, what 
we get is a result valid to 0(1/ d). Even in this case the 1/d method gives reasonable 
results that compare very favorably with other methods (see Ref. || and the references 
cited therein). 



3 



The free energy is now obtained by the standard trace over the fluctuating fields as 
well as the auxiliary fields and their Lagrangian multipliers 

T = -kT In y ... Jvu(z)Vp(z)Vg(z,z')VB(z,z')e~ m . (4) 

Once we have an explicit form for the free energy we can get an equation of state for the 
self-interacting chain, connecting the external stretching force acting on the chain with 
its elongation, from 

-^ = ((u z (L)-u z (0))) = (L-L ), (5) 

where L is the length of the chain after and Lq before the onset of the external stretching 
force. 



4 The Mean Field Solution 

Before integrating over all the fluctuating fields let us investigate the mean-field solutions 
of the mesoscopic Hamiltonian Eq. |3]. Let us first imagine we have no external tension 
applied to the chain, i.e. / = 0. The presence of the intersegment interactions, however, 
acts as an effective stretching force by itself. Let us see how that happens. The mean-field 
solution for this case is obtained by minimizing the Hamiltonian Eq. |3| and assuming that 
all the fields are constant: 

u zz = C, and p(z) = 0. (6) 

Thus since 

Hq = \\ J dz( 2 + \Jj dzdz'V(B(z,z')) + 

+ \ J J dzdz'g(z, z') (B(z, z') -(z- z' f (1 + C) 2 ) , (7) 
by minimizing with respect to ( we obtain the equation of state in the form 

where we introduced 

5X = J dz g(z, z ) (z — z') 2 . (9) 
The mean field minimization with respect to the auxiliary fields gives 

B(z,z>) = ( z - z f(l + (f 

g(z,z') = -d B V(B(z,z')). (10) 

We note at this point that the above mean field equations are highly and essentially 
non-linear. First of all g(z,z') is a non-linear functional of B(z,z r ), Eq. 0, and SX is 
determined from a solution to Eq. ^. 
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Thus on this level we see that the intersegment interactions stretch the chain in a way 
similar to an external force, leading to what one could call a Hooke's law of the form Eq. 
[S[ If we now add a real external tension (/) to the chain the mean-field ansatz would 
assume the form 

u zz = C + 5((f), and p(z) = 0. (11) 
The corresponding Hamiltonian in this case is 

Ho = |A / dz(( + 5() 2 -f J dz6( + \ J J dzdz'V(B(z, z')) + 

+ \ J J dzdz'g{ Zl z') (B(z, z') -{z- z'f (1 + C + Kf) • (12) 

Obviously we have coupled the stretching force only to the deformation (5() after the 
intrinsic deformation (() set by the intersegment interactions has been already established. 
Minimizing with respect to 5( we now get 

*f = A (13) 

and minimization with respect to the other variables gives 

B(z,z') = (z-z') 2 (l + C + 5C) 2 

g(z,z') = -d B V(B(z,z')). (14) 
Putting the two results together, the mean-field theory thus gives for the total deformation 

Again it is quite obvious that the intersegment interactions act in a way similar to an 
additional stretching force. Since the intersegment interactions make an additive contri- 
bution to /, see Eq. [TJ|, they obviously just displace the mean field minimum around 
which the system fluctuates. 



5 Fluctuations 

Now that we have the mean field solution for the case with external stretching force as 
well as intersegment interactions, we can expand the mesoscopic Hamiltonian around the 
mean field and evaluate also the effect of thermal fluctuations. To second order this 
expansion yields 

n = 




dzdz'g(z, z') ((u z (z) - u z (z')f + (p(z) - p(z')f) 



(16) 
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where Tio is given by Eq. OL In order to evaluate the functional integral corresponding 
to this effective Hamiltonian we first of all develop u z (z) and p(z) in the last term of 
the above equation into a Taylor series with an argument z — z'. This means that all the 
properties of the chain are homogeneous and depend only on z — z' M. We are thus trying 
to account for the longest length scale effects of the interaction terms on the properties of 
the semiflexible chain (an equivalent procedure would be to look at the lowest wave vector 
dependence of the Hamiltonian in the Fourier space). Both devices are consistent with a 
macroscopic character of the approach advocated here. The result of this expansion is as 
follows: 



\ \ \ dzdz'g(z, z') (u z {z) - u z (z')f = ±5\ / dz ( ] + . . . , (17) 



'du z (zY 

dz 



and 

1 f Uzdz ! g{z,z'){p{z)-p{z')f = \5\Jdz[^-^ - l 5Kc J dz [^^-\ +... 

(18) 

where the dots stand for higher derivative terms and we introduced 

5K C = ±( dz'g(z,z')(z-z') 4 . (19) 



'dp{z)\ 2 lST s f,(d 2 p{z) 



The complete Hamiltonian, including the mean field part as well as the contribution of 
fluctuations around the mean field now becomes 



H = \\Jdz(C + Kf~f Jdz5C, + \J Jdzdz'V{B{z,z')) + 



'd 2 p(z) 
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+ iA«/^j +^*H^J ■ (20) 

where we introduced the following renormalized elastic constants and a rescaled stretching 
force 

\W = X-5X 
K { C R) = Kc + SKc 

f = \(C + 5C)-5\ = 1— ^. (21) 

1 A 

This functional integral can be evaluated exactly for the harmonic variables u z (z) and p(z) 
assuming that we can ignore the end effects. The evaluation of the functional integral 
over non-harmonic degress of freedom , i.e. for the two auxiliary fields B(z — z') and 
g{z — z'), is dealt with on the saddle point level which constitutes the 1/d approximation 



(for details see Ref. 0) and leads to Eqs. (L4 
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The free energy of the chain can therefore be obtained in the form 



T = - kT In J ... J Vu(z)Vp(z)Vg(z, z')VB{z, z')e~ m = 
= \\ J dz(( + 5(f + f J dz5( + \J J dzdz'V(B(z, z')) + 

+ f lndet^)^) + fcT indet^^-/^). (22) 

The fluctuation determinants can be evaluated in the Fourier space by the standard 
methods ||. Since we have derived an explicit form for the free energy we can thus 
obtain the equation of state from Eq. in the form 

P- — - kT X J— 

Obviously the second term on the r.h.s. of the above equation comes from the transverse 
(p(z)) fluctuations and is thus entropic in origin while the last term is the mean field 
stretching term. We could also call them entropic and enthalpic stretching terms. 

We see immediately that in the case of no intersegment interactions or if the range 
of these interactions goes to zero (both of these cases leading to 5\ = 0, see below) the 
above equation of state reduces exactly to the one obtained by Odijk [JD]] and Ha and 
Thirumalai [J]. A similar equation of state has also been obtained by Marko and Siggia 
T0[ for the case of a chain with intersegment interactions except that the stretching part 



was added in by hand and that the bending and stretching moduli renormalization were 
not coupled as they are in Eq. |23|. 



6 Elastic Moduli Renormalization 

We now assume that in a uni-univalent salt solution the intersegment interaction potential 
is purely repulsive and of a screened Debye-Huckel form, i.e. 



v(r(^))= k ^ ^ { rf; {z) ;f h , 

a z \r(z) — r z 



where Ib is the Bjerrum length, a is the effective separation between the charges along the 
chain and k is the inverse Debye length. With this intersegment potential and assuming 
the mean field form for B{z) Eq. |14] one gets for the interaction driven changes in the 
stretching and bending moduli the following relations 



SK C = -- 2 



ds s 2 V'(B(s)) = - Ei(-m)) 

A 6 a 2 
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where Ei(x) is the standard integral exponent function and we introduced the local 
stretching parameter A as 

B(z,z>) ((rW-r(/)) 2 ) 2 (X + fV 

Obviously the renormalizations of the elastic moduli depend on the magnitude of the 
intersegment interactions (described by a) as well as on their range (set by the Debye 
length k^ 1 ). 

If the chain is inextensible, A — > oo, then A — > 1 , the renormalization of the 
stretching modulus becomes irrelevant and the renormalization of the bending modulus 
(second equation in Eqs. |25|) becomes exactly the Odijk-Skolnick-Fixman result @], as in 
deed it should. Also one realizes that experimentally A 3> / and thus one usually has 

A ~ A/(A - 5X) > 1. (27) 

The above relation Eq. ^ can be thus viewed as a generalization of the OSF result 
for the bending as well as stretching moduli. Since in the presence of the intersegment 
interactions A is a function of SX we have a very non-linear system of equations to solve. 
The solution would give us simultaneously the renormalization of the bending as well as 
the stretching moduli. 

Without even solving this set of equations we already know that the intersegment re- 
pulsions renormalize the bending and stretching moduli in the opposite directions, see Eq. 
^lj While the bending modulus increases, the stretching modulus decreases. The simple 
physical reasons for this were already outlined in the Introduction. This resolves com- 
pletely the conundrum observed in experimental studies of DNA stretching and bending 
elasticity [§]. 

One can furthermore examine the relation between renormalized, in effect measured, 
bending and stretching moduli in more detail. For their bare values we should have from 
the standard macroscopic theory of elasticity |J the result 

K c = \XR\ (28) 

where R is the radius of the molecule (for DNA R is between 10 A and 4 A , corresponding 
to phosphate and major groove radii). Obviously for the renormalized values Kq and 
this relation does not hold anymore. Instead we obtain the following relation between 
renormalized stretching and renormalized bending moduli 

= lA^i? 2 + (i + 2M) 2 A(1 - Ei{-Ka)j) . (29) 

The relation Eq.^ is thus valid only asymptotically as the range and/or the magnitude of 
interactions becomes very small. Any polyelectrolyte in the range of conditions where Eq. 



24| is valid should thus behave as a classical macroscopic cylinder if we take the coupled 
bending-stretching moduli renormalization due to the intersegment interactions properly 
into account. 
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7 Comparison with Experiment 



We can now fit the expressions Eq.|25] to the recent experiments by Baumann et al. f| 
where they measure the simultaneous dependence of the bending as well as stretching 
moduli of DNA on the ionic strength in a uni-uni valent electrolyte. The renormalized 
moduli for different values of the added uni-uni valent electrolyte are extracted from the 
fit of the experimental entropic and stretched regimes to the equation of state Eq. |23|, 
just below the overstretching transition. 

Fitting the dependence of the renormalized bending and stretching moduli on the 
inverse Debye screening length to Eqs. ^5], we can obtain both the bare bending and 
stretching moduli as well as the inverse line charge density a. Unbiased fits give for 
the bare elastic modulus 2.1 x 10 4 pN A (corresponding to lp = 511 A) and the bare 
stretching modulus A = 1511 pN. The values of the two bare moduli are completely 
consistent with Eq. [28| considering the fact that for DNA the radius R lies somewhere 
between the outer phosphate radius (~ 10 A) and the inner radius of the grooves (~ 4 A). 

From the fit to Eq.^5| we also obtain a consistent estimate (in the sense that it should 
fit the stretching as well as the bending modulus data) for a ~ 2 A. This estimate is not 
particularly accurate because of the large scatter present in the data, see Fig. ??, and 
because of the number of the fitting parameters. The experimental scatter is probably 
due to the fact that the regime between entropic and enthalpic stretching is quite narrow 
and a reliable estimate for the stretching modulus which can only come from this regime 
is thus difficult to obtain. Experiments are currently under way to gather a much more 
accurate set of data for the two elastic moduli [12[ . 



8 Discussion 

The theory presented above giving the coupling between stretching and bending moduli 
renormalization in the presence of finite range intersegment interactions seems to work 
reasonably when compared to experiments. There is of course no a priori reason for this 
and considering all the approximations that underpin the main results, Eqs. |2^, is perhaps 
quite surprising. 

First of all the result Eq. ^5] is formaly valid only in the limit of either a very stiff 
chain or very large external tension. If the stiffness or the tension are finite, we know 
|| that the extended configuration of the chain lying at the bottom of the Monge-like 
parameterization is unstable against thermal fluctuations. Nevertheless the OSF limit 
apepars to be stable [^] even in the regime of vanishing stretching and we expect (without 
any proof at this point) that the renormalizations Eqs. |25] will remain likewise. It would 
neverteless be appropriate to derive a more sophisticated theory, somewhat along the 
lines of but without the unstretchability constraint, that would be able to describe 
the equation of state for a stretchable self interacting semiflexible polymer chain for the 
whole range of stretching forces. 
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When calculating the fluctuation contribution to the free energy we assumed that 
B(z,z') and g(z,z') are still given by their mean field expressions, Eqs. H. There is in 
deed a fluctuation contribution to the auxiliary fields but it is in general small and would 
not fundamentally change the results derived above. It would however make the numerics 
more cumbersome. 

The numbers extracted for the effective charge density along DNA bear no resemblance 
to the Manning theory where the value appropriate for the effective spacing between the 
charges is a ~ 7 A. However the fits to the experimental data are not particulary stiff 
and we could certainly push the extracted numbers towards this value if we wanted. But 
considering the large scatter in the experimental data especially for the stretching modulus 
we produced an unbiased fit and intend to refine it as better data become available. 

Also we note that the form of the equation of state Eq. ^3| coincides with the one used 
to extract the values for the elastic constants || only in the no-intersegment-interactions 
limit. In general it gives corrections to this limit dependent on the salt concentration. 
These corrections should mostly affect the bending modulus which is obtained from the 
fit to the v7 dependent (entropic) part of the equation of state, and only marginally the 
stretching modulus, which is obtained from the fit to the linear part of the equation of 
state. This is another source that could eventually change the fitted parameters. 

Apart from all these numerical shortcomings and problems we advocate a fundamental 
change in the perspective regarding the dependence of the elastic moduli of DNA on the 
ionic conditions. Not only are the persistence length as well as the stretching modulus 
dependent on the salt conditions in the solution, they are fundamentally coupled. This 
is a consequence of the fact that as soon as the intersegment interaction potential is of 
finite range the stretching and bending themselves become coupled. This is most clearly 
exemplified by the exact form of the equation of state Eq. |23|. We believe that future 
work on the elasticity of DNA and similar (bio)polymers will have to take this fact into 
account. 
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10 Figure captions: 



Fig. 1 The physics of elastic moduli renormalization in the presence of finite range 
intersegment interactions. Repulsive intersegment interactions make the chain more dif- 
ficult to bend because of deminished effective separation between neighboring segments 
{V < L). However they also make stretching easier because they incerase the average 
length of the segments and thus diminish the curvature of the interaction energy of the 
segments (symbolically depicted as beads). 

Fig. 2 A highly stretched polymer chain. The average direction of the chain is along z 
axis, which corresponds also to the stretching axis and is set by the direction of external 
force / acting at both ends, and the bending deformation is perpendicular to it. 

Fig. 3 Experimental points taken from Table 1. of Baumann et al. and fits using 
the functional form of the ionic strength dependence from Eq.|25|. o are the measurements 
of persistence length ( left scale), • are the measurements of the stretch modulus (right 
scale) for ionic strengths between 1.86 and 586 mM. The functional form of the two moduli 
seems to fit the data quite well. 
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